Low-energy behavior of spin-liquid electron spectral functions 
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We calculate the electron spectral function for a spin-liquid with a spinon Fermi surface and a 
Dirac spin-liquid. Calculations are based upon the slave-rotor mean-field theory. We consider the 
effect of gauge fluctuations using a simple model and find the behavior is not strongly modified. 
The results, distinct from conventional Mott insulator or band theory predictions, suggest that 
measuring the spectral function e.g. via ARPES could help in the experimental verification and 
characterization of spin liquids. 



I. INTRODUCTION 

Anderson first proposed the resonating valence bond 
state in 1973 - an example of a spin liquid: states 
with no preferred spin orientation and no broken sym- 
metries. Such states may arise in frustrated spin systems 
such as the Heisenberg antiferromagnet on a kagome lat- 
tice. Only more recently have material candidates with 
promising experimental indications been found. The or- 
ganic material. k-(BEDT-TTF)2Cu2(CN)3 and a related 
material based on EtMe3Sb[Pd(dmit)2]2 molecules^ ex- 
hibit spin-liquid properties in the Mott insulating state, 
i.e. a lack of spin-order even down to very low temper- 
atures. Further, their linear specific heat behaviour at 
low temperatures points toward the presence of a finite 
density of states of gapless excitations, whereas thermal 
transport measurements on the latter compound suggest 
these excitations are mobile^. Another compound Her- 
bertsmithite ZnCu3(OH)gCl2 appears to have spin-liquid 
properties'^— with gapless spin excitations. 

The organic materials are "weak" Mott insulators 
which can be driven into a metallic phase under pres- 
sure. It has been proposed that the gapless excitations 
in such a weak Mott insulator can be modeled in terms 
of a Fermi sea of spinons - fermions that carry the spin, 
but not the charge, of the electron^ii. Indeed, recent 
work has established strong evidence that such a spinon- 
Fermi sea state occurs in a triangular lattice Heisenberg 
model with four-site ring exchange, a model appropri- 
ate for the organic Mott insulators^iS. An approach em- 
ploying Fermionic spinons has also been proposed as an 
explanation for the gapless behavior observed in the pu- 
tative spin liquid phase in Herbertsmithite, except with 
the spinons filling a Dirac sesi^—. 

In addition to thermodynamic and transport exper- 
iments, signatures for such spin liquids might be ex- 
pected in single electron properties such as those probed 
by angle-resolved photoemission spectroscopy (ARPES) 
or electron tunneling. In a conventional Mott insulator 
which exhibits antiferromagnetic order, the unit cell is 
doubled and the single hole spectrum can be described in 
terms of conventional band theory. For example, in high 
temperature superconductors, the holes occupy a band 
with a minimum at (7r/2, 7r/2)^ii^. For a spin liquid the 
unit cell is not doubled, and band theory predicts a metal 



with a half-filled band, which completely fails to describes 
the Mott insulator. It follows that the single hole spec- 
trum cannot have a conventional band-like description. 
Instead, within a slave-boson mean-field theory, the low- 
lying physical electron excitations are composites of a 
boson and the fermionic spinons and hence the spectral 
function behavior depends on the spinon dispersion. 

In this paper we work out the single electron spec- 
tral function for two models with different spinon spec- 
tra. Our goal is to stimulate experimental efforts to per- 
form ARPES on spin-liquid candidate materials. The 
first model is one where the spinons form a Fermi surface 
(called a SFS spin liquid). There is considerable evidence 
that this state may describe the organic spin liquid ma- 
terials. A second model is the Dirac spin liquid (DSL) 
state. In particular the tt-Aux phase on the kagome lat- 
tice has been proposed to describe HerbertsmithitoiSrJ^. 
While more recent work supports a gapped spin liquid 
as the ground state of the S = 1/2 antiferromagnetic 
Heisenberg model on the kagome lattice^^, it is believed 
that many states very close in energy are competitive and 
the DSL remains an interesting state. 



II. SPECTRAL FUNCTION BEHAVIOR 



The Mott transition between a Fermi-liquid metal and 
a spinon Fermi sea spin-liquid insulator can be described 
with a slave-rotor formalisnii^. Here the electron cre- 
ation operator is written as a product of a charged spin- 
less boson hi ("chargon") and a spin-1/2 charge-neutral 
"spinon" fi^,, i.e. fi^h. 



At the mean-field level, in the metallic phase the 
bosons condense and the spinon Fermi surface becomes 
an electron Fermi surface with a quasiparticle weight pro- 
portional to I (6,;) p. In the Mott insulator, the bosons 
are gapped while the spinons have a free hopping Hamil- 
tonian. An emergent gauge field = (ao,a) cou- 
ples to both spinon and boson fields fa and 6, and 
the resulting theoryii can be described by the action 
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FIG. 1: (Color online). Spectral function for a SFS spin- 
liquid at constant energy (satisfying \uj\ > + A): excited 
electrons appear as a ring with radius Kf in the Brillouin 
zone. The ring broadens as a function of energy — ring width 
goes as 2(5fcmax ~ 1\j2m\\bj + Aj. Maximum intensity is at 
|k| = Kf where ylo(KF,i.j) = ^2m'\i^ + A|. 



S = dr J d'^rC, where 
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Cb + Cf + Cg, 

|(9, + ^ao)6p + IKV + ^a)6|2+A2|6|^ 
2m/ 



(1) 



Here is the low-energy boson velocity, A the insulating 
gap and m f the effective spinon mass. Beyond mean- field 
the gauge-field fluctuations have to be considered, which 
will be discussed at the end of this paper. 

At the mean- field level, the electron correlator factor- 
izes locally (in space and time) into boson and spinon 
correlators. The retarded Green's function for the elec- 
tron is then given by a convolution of boson and spinon 
Green's functions. At zero temperature we getii 
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where Ck-q 



is the dispersion of the a spinon band, flq 
is the boson dispersion and / is the Fermi distribution 
function. 

The form factor (k — q, q) depends on the spinon 
and boson wavefunctions V'" a-nd 0, 



Af„(k-q, q) 



iJCiy^ - q)VC(k - q)0™(q)C(q), 



where m runs over each atom in the basis. The derivation 
and evaluation of this expression is given in Appendix |B] 
— at low-energies this form factor is simply a constant 
and does not change the qualitative result. 

We obtain the electron spectral function as usual us- 
ing Ao(k, w) = — ^IiyiGq (k, Lu) . For electron excitations 
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FIG. 2: Spectral function for a SFS spin-liquid at fixed mo- 
mentum k = Kf + (5k where 5k is parallel to Kf : we expect 
a square root dependence when varying lj (past the threshold 
energyof |a;o| = |^+A). 



(negative uj) we keep only the second term in Eq. [21 
At low-energies, the universal properties of the spectral 
function are determined only by the spinon and boson 
dispersions. For simplicity the boson dispersion is as- 
sumed quadratic. 
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(3) 



where A is the insulating gap and mf, = A/?;^. Here 
(5q is the boson momentum q measured from the band 
minimum or minima, which are located at Qb. 

We consider two cases for the spinon dispersion. In the 
first there is a spinon Fermi surface and the second is a 
Dirac dispersion with a vanishing density of states at the 
chemical potential. 



A. Spin liquid with spinon Fermi surface (SFS) 

In the model proposed for the organic spin-liquid^ii, 
the spinous have a Fermi surface and the boson minimum 
is located at Qb = 0. At low energies we obtain 



Ao{Kf +dk,io) 




2mh 



(4) 



where Sk is along Kp and = — —. This has a 

square-root dependence on \iu\ above the threshold exci- 
tation energy, jwo] = ^ — h A. Spectral function mea- 
surements could detect such behavior, say using ARPES. 
When examining electronic excitations at constant en- 
ergy (satisfying > ^--|-A), one would measure a ring 
with radius Kf in the Brillouin zone, see Fig. [TJ The ring 
broadens as a function of energy, with ring width going 
as 2(5fcniax = 2\j2rrn,\jj + A|. The maximum intensity is 
at |k| = Kf where Ao(KF,a;) ^2m!\bj + A|. 

If the experiment is done instead at fixed momentum 
k, we expect a square root dependence when varying to 
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FIG. 3: (Color online). Spectral function for a DSL at con- 
stant energy (satisfying jcj| > + A): excitations are ex- 
pected at certain points {Ke} in the Brillouin zone. The 
dot radius size is a function of energy 5fcmax ~ \j2m\y\u) + A| 
and maximum intensity is at the center where j4o(Ke,Cij) = 
2An^ 1^ + ^1 ■ Note that here a kagome lattice with the 7r-flux 
state is shown for concreteness; however, a generic DSL will 
have the same behavior for low-energy excitations. 



once electrons have enough energy to overcome the gap 
A, see Fig. [H 
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FIG. 4: Spectral function for a DSL at fixed momentum k: we 
expect a linear dependence when varying to past the threshold 
energy of |a;o| = ^ + A with slope -gf^- 



with energy Jfcmax = ■\/2mf,|a; -I- A| and maximum in- 
tensity would be observed at the center of these dots, 
Ao(Ke,w) = 2^J^I^ + ^1- If ^ measurement is done at 
fixed momentum k, we expect instead a linear depen- 
dence when varying to with slope — ^^^.^ , see Fig. U) 

III. TT-FLUX STATE ON THE KAGOME 
LATTICE 



B. Dirac spin liquid (DSL) 

For spinous with a Dirac spectrum, the dispersion is 
linear to lowest order in momentum 



(5) 



where v is the spinon velocity at the Dirac point and (5k 
is k measured from the Dirac point. 

After integration (details in Appendix the spectral 
function becomes 



Ao(Ke + (5k. 
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The constant M' is the form factor contribution (defined 
and calculated in Appendices \^ and [B| which simplifies 
considerably near the Dirac point. 

The result has a linear dependence on |a;| past the 
threshold excitation energy, which we expect to see 
around a set of points of electron excitations Ke in the 
Brillouin zone. For concreteness we show this behavior 
for a TT-flux state on the kagome lattice (details regarding 
this model arc in the next section). 

In the hexagonal Brillouin zone of the kagome lattice, 
electron excitations are expected from eight points, see 
Fig. |3l Here we plot the spectral function at fixed en- 
ergy. Once there is enough energy to overcome the gap 
A, excitations appear at these points and broaden as a 
function of energy. The radius of these dots increases 



As the results in the previous section are fairly gen- 
eral, we do not expect them to be strongly modified by 
microscopic details or the model used. The form-factor 
Ma{\i — q, q) however, can be complicated and depends 
explicitly on band structure and lattice details. Here we 
pick a particular model to do an analytical expansion and 
calculation to make sure the form factor does not vanish 
or show additional singularities. 

Further, DSL models depend on a gauge choice that 
initially appears to break translation invariance. While 
spinon parameters do depend on the gauge chosen, the 
spectral function as a convolution of boson and spinon 
momenta is gauge-invariant (as expected for a physical 
observable). In our model we see explicitly that neither 
the gauge choice nor microscopic details affect our results. 

The TT-flux state on the kagome lattice was proposed 
as a candidate state for the Dirac spin liquid.— In this 
model the tight-binding Hamiltonian describes a tt-Aux 
through every hexagon and zero flux through each trian- 
gle on the kagome lattice (see Fig. [5]). As this doubles 
the unit cell, the spinon and boson dispersions have six 
bands. With nearest-neighbor hopping t, we obtain the 
following dispersion (in units where the magnitude of the 
lattice spacing is set to 1/2) 
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2t (doubly degenerate) 
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At half-filling, the spinon chemical potential [i = 
(-\/3 — \)t is located between the third Ej^- and fourth 
i?4- + bands which touch at two Dirac nodes k = 
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FIG. 5: (Color online). Kagome lattice with a doubled unit 
cell (translation vectors ri and r2). Dashed lines correspond 
to bonds with negative hopping while unbroken lines corre- 
spond to bonds with positive hopping — such that hexagons 
enclose a vr-flux while triangles have zero flux. The band 
structure given by this hopping model (plotted along the line 
kx = 0) shows six bands (the top band is twofold degenerate). 
Spinon excitations are from the Dirac points at half-filling 
fj, = (V^ — l)t while boson excitations are from the lowest 
band minima. 



(0, ±7r/A/3)(Fig. O. For electron excitations, the spinon 
band a is just the fourth band. Bosons excitations are 
from the minima of the lowest band E_ .^_ — points 
k = {±(7r/3,0),±(27r/3,7r/\/3)} in the Brillouin zone. 

For low-energy excitations we expand to second order 
around a band minimum for bosons and a Dirac point 
for spinous, 



A = 



4V6 



(7) 



.lL\Sk\ + ^idkr + 0{Sk') 



where 5q and (5k are momenta measured from the boson 
minimum and Dirac point respectively. Here the nearest- 
neighbor hopping t has been replaced hy tf and tb for 
effective spinon and boson hoppings obtained from the 
self-consistent mean-field solution. So in this model the 
boson mass and spinon velocity are mi, = 2\/6/tb and 
V ^ tf/V2. 

Since low energy excitations are composed of spinous 
near the Fermi energy and bosons near the band bottom, 
electron momenta will be the sum of both spinon and 
boson momenta. All possible combinations give rise to 
excitations centered around eight points in our halved 
Brillouin zone 



K, 



= (±- ± — 



27r \ /27r 27r , 



(The last two points can be translated to four other 
equivalent Kg points in the larger Brillouin zone using 
the original reciprocal lattice vectors.) 

Given the equivalent behavior of excitations in this sys- 
tem, it suffices to calculate the form-factor around one 
particular excitation where we find M' w 0.6 (see Ap- 
pendix |B|) . 
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FIG. 6: (Color online). A comparison in our DSL vr-flux 
model of the general linear prediction from Eq. [S] (straight 
line) with the actual mean-field behavior given by Eq. lAll 
(dotted lines). Here we choose low values of momenta close 
to one excitation point k ~ (vr/S, 7r/\/3). With tb — tf — 1, 
we compare Ae(k,ijj) for two values of cSk = (0.1,0.1) and 
(5k — (0.2,0.2). We see the analytical expression is a good 
approximation to the mean-field expression at low energies, 
and is better for smaller values of (5k as expected. 



We would like to compare the general linear predic- 
tion from Eq. [6] with the actual mean-field behavior 
of the spectral function in this DSL model. Here we 
choose low values of momenta close to one excitation 
point, k (tt/S, 7r/-\/3). Setting tb = tf = 1, we 
compare Ae(k, w) for two values of Sk = (0.1,0.1) and 
Sk = (0.2,0.2) (see Fig. [S]). We see the analytical ex- 
pression is a good approximation to the mean-field ex- 
pression at low energies, and is better for smaller values 
of Sk as expected. 



IV. EFFECT OF FLUCTUATIONS 

Beyond mean-field, the dominant effect of gauge-field 
fluctuations is to bind the spinon and boson into an elec- 
tron. Here we consider a simple model for these fluctua- 
tions to find the result does not change qualitatively for 
weak interaction strengths. 

It is possible to approximate the attraction with an 
effective short-range interactio n-'^^'^^ of the form 



U>0. 



This changes the electron Green's function to, 
Ge{k,uj) ^ 



(9) 



where Go(k,io) is the mean-field expression given in Eq. 
m We denote its real and imaginary parts G'o(k, w) = 
G'o{Koj)+iG'^{k,io). 

Having already calculated the imaginary part in Sec- 
tionHH we need only to further calculate Go(k, uj). In the 
following we focus on the Dirac spin-liquid. As detailed 
in Appendix [Cl in this case one finds. 



Uc \ qc \ M 



(10) 
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r > 2 

Here we have defined w = (w + |^ + A)/v and also 
Uc = M'gc/2A, with Qc a high momentum cutoff. 

The electron spectral function in the presence of the 
attractive interaction between the spinous and chargons 
can now be obtained by taking the imaginary part of 
Ge(k,w) in Eq. 1 



ImGe(k;, lu) 



G" 



(1 



Notice that there is a pole only when Gq = (i.e. when 
(D > 0) and Gq = so a transition occurs at the 

binding energy Uc- 

When U < Uc, there is no bound state, so the pres- 
ence of the spinon-chargon interaction merely alters the 
original linear behavior of ^o(k, w). To leading order in 
uj the interacting electron spectral function becomes, 



Ae(k,a;) 



Ao{k,io) 



^o(k,a;) 

il-U/Ucr' 

i.e. Ae{\i,Lu) is still linear with an enhanced slope. 

When U > Uc, there exists a bound state between the 
spinon and chargon (forming an electron) for Gg = l/U, 
implying an electron state at energy uj satisfying, 



)ln 



= qc[i- 



U 



(11) 



The electron spectral function has a delta-function at 
this value of uj and a small incoherent part where a) < 0. 
The above equation can be solved numerically, e.g. for 
U = 2Uc and Qc = 1 we get Co = 0.4. In this case uj = 
ujv — — A, i.e. the (negative) energy of the bound 
state within the gap decreases with Sk^ /2mb and has its 
minimum on the Fermi surface. 

Hence, for weak interaction strength U, Ae(k.,uj) has 
the same dependence on uj and Sk. as the mean-field result 
— only the coefficient is modified. 



V. DISCUSSION 

We present particular signatures in the electronic spec- 
tral function of spin-liquids that would distinguish such 
materials from a conventional Mott insulator or a con- 
ventional band theory description. Using slave-rotor the- 
ory at the mean-field level, we consider two spin-liquid 
models: one with a spinon Fermi surface and another 
with a Dirac spinon sea, and extract the associated elec- 
tron spectral functions. We also explore the role of 
gauge fluctuations employing a simple model consisting 
of an attractive spinon-chargon interaction. Solving this 
simple model indicates that weak interactions do not 
change the qualitative behavior of the mean-field spec- 
tral functions, but above a threshold interaction strength 



a spinon-chargon bound state (an electron) forms. The 
energy of the bound state peels off from the continuum 
energy threshold, leading to a delta-function in the elec- 
tron spectral function inside the gap. We hope this work 
will encourage experimental efforts to perform ARPES 
on spin-liquid candidate materials. 

After completion of this paper, we learned that Pujari 
and Lawlcr also studied the ARPES signature of the DSL 
state on the kagome lattice, to find the linear onset is 
unaffected by gauge fluctuations to first order—. 

We thank Andrew Potter and David Mross for help- 
ful discussions. This work was supported by the NSF 
under grants DMR-1104498 (P.A.L.) and DMR-1101912 
(M.P.A.F.) and by the Caltech Institute of Quantum In- 
formation and Matter, an NSF Physics Frontiers Center 
with support of the Gordon and Betty Moore Foundation 
(M.P.A.F.). 



Appendix A: DSL spectral function integral 

Here we evaluate the spectral function 
Ao{k,io) = ^-i-(Af„(k-q,q)^(c.-e£_q-f^q) 



+M„(k-q,q)(5(^-e£_q + l^q)), (Al) 

for a Dirac spin liquid. 

As the boson wavcfunctions are quadratic at the band 
minima, A/(k — q, q) depends only on spinon momenta. 
Taking a to be the spinon band below the Fermi energy, 

Ao(k,uj) 



2^/6 2mb 



where we changed coordinates to q' ~ Sk — (5q, momenta 
about the Dirac points. 

As the low-energy spinon wavefunctions do not depend 
linearly on |g'| but only on 9 — the direction of q' around 
the Dirac point — M{q') reduces to M{9) (independent 
of the magnitude of 6k — 6q). In polar coordinates the 
integral becomes 

Ao{k,Uj) = ^ l^\eM{e) dq'q'5{uj + v\q'\ 



vq'^ {5k - q'f 



2^/6 2mf, 

c'2-K pqc 



A) 



1 r^'" /""^^ Sk^ 



lo Jo 



2At;2 ' 2mb 



f)k'^ 

\u; + —\>A, (A2) 
Znib 
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where we define the constant M' = d9M (6) 



Appendix B: Form-factor contribution M' 

MQ(k— q, q), the "form factor" given in Section HIl 
depends on the spinon and boson wavefunctions tp" and 
(j) (which have m indices for each atom in the basis) — 
so is an m X m matrix, 



MQ^,„„(k- q, q) = 



^ gik-(r„-r„) 



i=0,3 



Since the tt-Aux state has a doubled unit cell, we sum over 
Z = 0, 3 to double the number of atoms as m,n = 1, 2, 3 
in the original basis. 

Resolution on an atomic level is uneccssary so we take 
the trace, 



M„(k-q, q) 



TrA4,,„„(k-q, q) 
^ Ma^mm (k - q,q) 




To simplify the expression further, we use certain prop- 
erties of the model that also demonstrate the explicit 
translation invariance (in the original unit cell) of our 
model and hence its gauge invariance. 

Each electron excitation has contributions from the 
two spinon Dirac points. The Hamiltonian at the first 
Dirac point is the complex conjugate of the Hamiltonian 
at the second Dirac point, 

Hik = (0, tt/Vs)) ^H*{k= (0, -n/Vs)), 

Also, at low energies a spinon from the first Dirac 
point combines only with bosons from the two min- 
ima q « {(7r/3,0),(-27r/3,-7r/y3)}, while spinous 
from the second Dirac point combine only with 
bosons from a non-overlapping set of minima, q w 
{(— 7r/3, 0), (27r/3, 7r/-\/3)}. The boson wavefunctions 
that pair with a spinon from the first Dirac point is the 
complex conjugate of boson wavefunctions that pair with 
the second: 

((/),„+Kq)C(q))* = (/.„+Kq± (^,7r/%/3))C(q± (^r, V^/3)) 

These properties enable us to simplify the form fac- 
tor — as each excitation has two contributions that are 
complex conjugates, their imaginary part cancels — 

M„(k-q,q) 

= ^"n+i (k - q)VC (k - q)0™+i (q)0;„ (q) 

m.l 

= ^|^;^(k-q)n0™(q)p (Bl) 

771 

+Rc[ij^,+i(k - q)V'r(k - q)0m+i(q)</',*„(q)] 



FIG. 7: A plot of Ma(6) against 9, it is a smoothly varying 
function. 



As the form factor is now explicitly Hermitian, it pre- 
serves translation invariance in the original unit cell and 
also gauge-invariance. 

We now calculate M' for an excitation near Kg = 
(7r/3, 7r/-\/3) with a spinon from k- q « (0,7r/^/3). At 
low-energies, a is the fourth band and we expand the 
Hamiltonian to a subspace of the two bands near half- 
filhng. Computing M{d) (see Fig. [7|), this is 

M' 

deM{e) 

d0^A/„„„(0) 

d6 {\An{0)?\<t>m? 
m=l,2,3 

Rc[^3+m(^)C(^)03+mC]} 



27r 



+ Rc[^4(e)^i*(e)040t + V'5(e)V^2W050; +^6(e)^3(^)'/>603]} 

0.604 (B2) 

so the form factor M' suppresses the spectral function by 
about 27r/0.604 « 10. 

Hence, the form factor merely enhances or suppresses 
particular excitations as a smoothly varying function of 
angle around the Dirac point and is independent of lo or 
the magnitude of k in this expansion. 



Appendix C: Binding energy integral 

We can use a Kramers-Kronig relation to obtain the 
real part of our mean-field Green's function from Eq. [2] 



Go(k,a; 
1 

2U 



k-q 



CJ — ^■ 



k-q 
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Keeping only electron excitations and expanding at low For a; < 0: 
energies, we get 

S^J_-pf MjSk-Sq) \ f 1 \ 

(C2) _ r\^\-n dqq ^ dqq 

Changing again to q' = — (5q and using polar coordi- "^^ q + "'|cj|+77 1 + ^ 

nates this becomes i^, ~, f V \ /,-, ~, flc + i^ 

— \uj\ - ?/ - win [ — ] + qc- (\oj\ + ijj — ojin 

G'o(k,w) 

1 f r . . .^1 1 t V iw 



deM{e) j dq'q'r 



^1/ \ V 

q^ u\n I (C5) 



2AJ -7 V- + -k'|-0i + ^ + A^ 

f d0M{e) f dq'q'V I ) 

2 A J J \vq' + Lo + 2;;^ + A/ Combining both expressions, we obtain the real part of 

J, I / \ the Green's function 

^ ^ ' dmv — , (C3) 



where qc is a cutoff momenta far above the scale of low ^ 
energy excitations. ^ (^^(/^^ ^) = ( _ Cuhi ' + 

5k I Z/^"' 



To evaluate this integral, we look at a; = (a; + 



1->P^ (C6) 



A)/z; > and w < 0. = J_ M _ ii^ln I + ^ 

For w > 0: t/, 



q + LO J Jq \ q + ijj 

= qc- ciln ( ^^-^ ) (C4) where we define Uc = M'qc/2A. 
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